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1.0  INTRODUCTION 

This  report  documents  the  specific  dynamics  and  mobility  factors  within  a  U.S.  Army  ground  vehicle 
Design  of  Experiments  (DOE)  approach  to  designing  an  Occupant  Centric  Platform  (OOP)  using  Multi¬ 
disciplinary  Design  Optimization  (MDO).  The  objection  of  an  OOP  is  to  “ 

Understand  how  to  balance  vehicle  “protection”,  “performance”  &  “payload”  through  an 
integrated  survivability  approach  that  starts  with  occupant  protection. — [1] 

As  [2]  states,  the  problem  is  that  “We  design  vehicles  to  put  Soldiers  in  rather  than  designing  vehicles 
around  Soldiers.  Increasing  protection  levels  of  the  platforms  impacts  interior  volumes  reducing  mobility, 
maneuverability,  and  freedom  of  movement  for  occupants  and  leads  to  heavier  platforms.”  while  the 
challenge  is  to  “Formulate  a  S&T  program  to  make  improvements  to  existing  platforms  or  develop  new 
platforms  that  provide  appropriate  increased  protection  from  current  and  emerging  threats  and  optimal 
space  allocation  for  Soldiers  and  their  gear,  while  decreasing  platform  weight  and  maintaining  or 
increasing  maneuverability  during  full  spectrum  operations.”  This  encompasses  a  wide  array  of  subject 
areas  in  vehicle  design  that  includes  Survivability,  Mobility,  Thermal,  Power,  Dynamics,  Structures,  etc. 
As  explained  in  [3],  MDO  addresses  the  complex  interactions  between  all  of  these  areas. 

Multi-disciplinary  design  optimization  (MDO)  is  a  field  of  engineering  that  uses 
optimization  methods  to  solve  design  problems  incorporating  a  number  of  disciplines.  As 
defined  by  Prof.  Carlo  Poloni,  MDO  is  "the  art  of  finding  the  best  compromise".  It  is  also 
known  as  multidisciplinary  optimization  and  multidisciplinary  system  design  optimization 
(MSDO). 

MDO  allows  designers  to  incorporate  all  relevant  disciplines  simultaneously.  The 
optimum  of  the  simultaneous  problem  is  superior  to  the  design  found  by  optimizing  each 
discipline  sequentially,  since  it  can  exploit  the  interactions  between  the  disciplines. 
However,  including  all  disciplines  simultaneously  significantly  increases  the  complexity  of 
the  problem. 

These  techniques  have  been  used  in  a  number  of  fields,  including  automobile  design, 
naval  architecture,  electronics,  [[architecture],  computers,  and  electricity  distribution. 
However,  the  largest  number  of  applications  have  been  in  the  field  of  aerospace 
engineering,  such  as  aircraft  and  spacecraft  design.  For  example,  the  proposed  Boeing 
blended  wing  body  (BWB)  aircraft  concept  has  used  MDO  extensively  in  the  conceptual 
and  preliminary  design  stages.  The  disciplines  considered  in  the  BWB  design  are 
aerodynamics,  structural  analysis,  propulsion,  control  theory,  and  economics. — [3] 

A  wide  range  of  software  tools  are  unique  to  each  discipline.  Several  MDO  software  packages  are 
available,  see  [4]  and  [5]  as  examples,  that  provide  solutions  for  handling  and  incorporating  the  results 
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from  each.  As  stated  in  [4],  these  are  used  to  combine  cross-disciplinary  models  and  applications 
together  in  a  simulation  process  flow,  automate  their  execution,  explore  the  resulting  design  space,  and 
identify  the  optimal  design  parameters  subject  to  required  constraints.”  Within  this  framework  the  initial 
simulation  complexity  can  be  reduced  by  using  simplified  models.  As  optimized  results  are  received, 
refinements  can  be  made  to  adjust  resolution  as  required. 

The  objective  of  this  report  is  to  provide  simplified  dynamic  and  mobility  factors  that  reflect  pertinent 
vehicle  design  performance  in  these  areas  and  that  can  be  used  in  an  MDO  for  an  OCP  with  an 
underbody  blast  survivability  focus.  The  calculation  of  each  factor  is  described  below  and  based  on 
portable  scripts  that  fit  well  into  the  framework  described  above. 


2.0  DYNAMICS 

The  end  result  of  neglecting  safety  can  be  as  consequential  as  survivability  threats.  Resistance  to  vehicle 
rollover,  both  military  as  well  as  civilian,  has  to  be  considered  in  vehicle  design.  While  rollover  is  a 
dynamic  phenomenon,  simplified  measures  can  be  used  to  gauge  the  rollover  propensity.  One  of  these 
is  the  Static  Stability  Factor  (SSF).  As  noted  from  [6], 


“The  Static  Stability  Factor  (SSF)  of  a  vehicle  is  an  at-rest  calculation  of  its  rollover 
resistance  based  on  its  most  important  geometric  properties.  SSF  is  a  measure  of  how 
top-heavy  a  vehicle  is. 

A  vehicle's  SSF  is  calculated  using  the  formula  SSF=T/2H,  where  T  is  the  "track  width"  of 
the  vehicle  and  H  is  the  "height  of  the  center  of  gravity"  of  the  vehicle.  The  track  width  is 
the  distance  between  the  centers  of  the  right  and  left  tires  along  the  axle.  The  location  of 
the  center  of  gravity  is  measured  in  a  laboratory  to  determine  the  height  above  the 
ground  of  the  vehicle's  mass.  The  lower  the  SSF  number,  the  more  likely  the  vehicle  is  to 
roll  over  in  a  single-vehicle  crash.” 

A  description  of  the  specific  vehicle  parameters  described  above  are  shown  in  Figure  1. 


Figure  1  -  Static  Stability  Factors 


From  [6]  further  details  are  given  that  “The  SSF  metric  is  based  on  a  rigid-body  model  of  a  vehicle  sliding 
laterally  on  a  surface.  For  such  a  model,  the  point  of  incipient  rollover  occurs  when  the  sum  of  the  lateral 
forces  divided  by  the  weight  of  the  vehicle,  W,  is  greater  than  the  SSF:  sum  of  lateral  forces/W  >  SSF  (= 
T/2H).  For  a  vehicle  to  roll  over,  the  lateral  forces  must  be  sustained  for  a  sufficient  period  of  time.” 
Findings  from  [6]  also  indicate  that  the  SSF’s  relevance  to  rollover  is  “based  on  the  laws  of  physics  and 
captures  important  vehicle  characteristics  related  to  rollover”  while  “metrics  derived  from  dynamic  testing 
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are  needed  to  complement  static  measures”.  From  these  conclusions  SSF  is  determined  to  be  a 
satisfactory  simplified  model  to  predict  one  important  element  of  dynamics  that  relates  to  OCP  design 
parameters. 


2.1  SSF  CALCULATION 

A  Python  (see  [12])  script  was  developed  to  read  an  input  file  of  variable/value  pairs,  calculate  the  SSF, 
and  write  out  a  file  of  the  variable/value  pairs  including  the  SSF  factor  (see  Appendix  A.1  -  latw.py). 
Including  the  input  factors  in  the  output  file  provides  context  for  the  output  calculation  and  can  help 
diagnose  errors  during  initial  development.  Reading  an  input  file  and  writing  an  output  file  allows  the 
script  to  be  easily  adapted  to  any  MDO  software  package  described  in  section  1 .0. 

The  SSF  is  returned  from  a  call  to  the  dynutil  library  (see  Appendix  A.3  -  dynutil.py)  as 

ssf  =  dynutil . calc_SSF (v_t,  v_h) 


and  defined  as 


where  higher  values  represent  better  stability  (The  HMMWV  is  around  1.0).  Once  the  input  file  has  been 
created  the  script  is  called  as  “latw.py  [input  file]  [output  file]”.  An  example  is  shown 
below. 


latw.py  latin.w  latout.w 


2.2  SSF  INPUT  FILE 

The  format  of  the  input  file  is  simply  variable/value  pairs,  one  set  per  line,  separated  by  tabs  or  spaces. 
An  example  is  shown  below. 

latin.w _ 

#=============================================================================== 

# 

#  Description: 

#  Input  file  defining  variable  values  to  Calculate  Static  Stability 

#  Factor  (SSF) 

# 

#=============================================================================== 

#  t  #  track  width  [in.] 

#  h  #  CG  height  above  ground  [in.] 

I : - - ==: - == - 

I  use  tab  or  space  separators 

#=============================================================================== 

t  76 

h  36 


2.3  SSF  OUTPUT  FILE 

The  format  of  the  output  file  is  simply  variable/value  pairs,  one  set  per  line,  separated  by  tabs.  An 
example  is  shown  below. 
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latout.w 


h  36.000000 
t  76.000000 
SSF  1.055556 


3.0  MOBILITY 

There  is  no  more  direct  measure  of  mobility  than  the  determination  of  whether  a  vehicle  becomes 
immobile— it  cannot  move.  This  is  done  by  comparing  the  available  traction  to  total  resistance.  In  a 
cross-country  setting,  grades  (longitudinal  slope),  side  slope,  engine  power,  etc.  are  several  contributing 
factors  but  soil  strength  is  the  key  driver.  The  immobile  or  “NO-GO”  area  is  demonstrated  through  the 
use  of  the  mobility  map  and  speed  profile,  in  Figure  2,  where  the  speed  falls  below  the  average  speed  for 
less  than  100%  of  the  terrain  area.  Soil  specific  details  are  outlined  in  [8]  where  the  beginnings  of  soil- 
vehicle  models  are  discussed  along  with  several  metrics  for  directly  quantifying  the  ability  of  vehicles  to 
traverse  soft-soil  terrain.  One  of  these  is  Vehicle  Cone  Index  (VCI).  As  noted  from  [9], 


VCI  is  defined  as  the  minimum  soil  strength  necessary  for  a  self-propelled  vehicle  to 
consistently  make  a  prescribed  number  of  passes  in  track  without  becoming  immobilized. 
Historical  testing  usually  focused  on  measuring  the  minimum  soil  strength  required  for  a 
vehicle  to  make  one  pass  (VCI1 )  and/or  50  passes  (VCI50).  Using  VCI  measurements  for 
a  number  of  different  vehicles,  the  US  Army  ERDC  developed  a  set  of  empirical 
equations  that  predict  VCI1  and  VCI50  from  relatively  simple  vehicle  characteristics  (i.e., 
weight  and  running-gear  dimensions). 


A  description  of  the  specific  vehicle  parameters  that  influence  the  empirical  relationships  are  described  in 
Figure  3  [9], 


V 


Figure  2  -  Mobility  Overview 
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Nomenclature 

b 

average  tire  section  width  (inflated;  unloaded),  in. 

CF 

clearance  factor  (part  of  MI) 

CPF 

contact  pressure  factor  (part  of  MI) 

d 

average  tire  outside  diameter  (inflated;  unloaded),  in. 

DCF 

deflection  correction  factor 

EF 

engine  factor  (part  of  MI) 

GF 

grouser  factor  (part  of  MI) 

GYW 

gross  vehicle  weight,  lb 

h 

average  tire  section  height  (inflated;  unloaded),  in. 

hc 

vehicle  minimum  clearance  height,  in. 

m 

total  number  of  axles 

MI 

mobility  index 

MMPm 

newly  proposed  version  of  mean  maximum  pressure  (by  Maclaurin),  psi 

MMPr 

original  version  of  mean  maximum  pressure  (by  Rowland),  psi 

n 

average  number  of  tires  per  axle 

TEF 

traction  element  factor  (part  of  MI) 

TF 

transmission  factor  (part  of  MI) 

VCIi 

one-pass  vehicle  cone  index,  psi 

VCI50 

50-pass  vehicle  cone  index,  psi 

w 

average  axle  loading,  lb 

WF 

weight  factor  (part  of  MI) 

WLF 

wheel  load  factor  (part  of  MI) 

(5 

average  hard-surface  tire  deflection,  in. 

Figure  3  -  Vehicle  Cone  Index  Factors 


From  [10],  The  International  Society  for  Terrain-Vehicle  Systems  (ISTVS)  defines  VCI  as  one  of  its 
standard  terms.  Its  usage  is  also  described  in  [11],  From  these  conclusions  VCI  is  determined  to  be  a 
satisfactory  simplified  model  to  predict  one  important  element  of  mobility  that  relates  to  OCP  design 
parameters. 


3.1  VCI  CALCULATION 

A  Python  (see  [12])  script  was  developed  to  read  an  input  file  of  variable/value  pairs,  calculate  the  VCI, 
and  write  out  a  file  of  the  variable/value  pairs  including  the  VCI  factor  (see  Appendix  A. 2  -  vciw.py).  As 
sated  in  section  2.1 ,  including  the  input  factors  in  the  output  file  provides  context  for  the  output  calculation 
and  can  help  diagnose  errors  during  initial  development.  Reading  an  input  file  and  writing  an  output  file 
allows  the  script  to  be  easily  adapted  to  any  MDO  software  package  described  in  section  1 .0. 

The  VCI  is  returned  from  a  call  to  the  dynutil  library  (see  Appendix  A.3  -  dynutil.py)  as 

veil  =  dynutil . calc_VCIl (v_gvw,  v_m,  v_b ,  v_d,  v_h ,  v_hc ,  v_n ,  v_delt ,  v_CGF ,  v_CEF,  v_CTF) 


and  defined  as 
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Vehicle  Cone  Index  (one-pass,  wheeled  vehicle) 
VCIi  =  f(MI,  DCF) 


where 


f  (  39.2  \ 

I  MI  <  115  ->  VCIi  =  f  11.48  +  0.2MI  -  —  -  -  —j  DCF, 

(MI  >  115  ->  VCIi  =  (4.1MI°-446)DCF. 


Mobility  Index 

/  (CPF)(WF)  \ 

M,  =  (frEFKCR+WLF-CFj(EFMTF) 


Mobility  Index  Factors 


CPF  = 


w 

0.5ndb 


TEF  = 


10  +  b 
100 


WLF  = 


w 

2000 


GF  =  1  4-  0.05CGF,  where  CGF  =  1  if  tire  chains  are  used  or  0  if  not 

hp 

EF  =  1  +  0.05CEF,  where  CEF  =  1  if  PWR  <  10  —  or  0  if  not 

ton 

TF  =  1  +  0.05CXF,  where  CXF  =  1  if  manual  transmission  or  0  if  automatic 

w 

WF  =  CWF1  Yqqq  +  ^WF2, 

{w  2000  lbs.  *  Cwfi  =  0.553  and  C^F2  —  0, 

2000  <  w  <  13,500  lbs.  ->  CWF1  =  0.033  and  CWF2  =  1.050, 
13,500  ^  w  20,000  lbs.  —>  C^vfi  =  0.142  and  C^yp^  =  0.420, 
20,000  ^  w  31,500  lbs.  —>  C^vfi  =  0.278  and  C^yp^  =  3.115, 
31,500  ^  w  —  0.836  and  C^yF2  —  20.686. 


Deflection  Correction  Factor 


DCF  = 


0.25 


where  lower  values  represent  better  mobility  (The  original  HMMWV  is  around  25.0).  Once  the  input  file 
has  been  created  the  script  is  called  as  “vciw.py  [input  file]  [output  file]”.  An  example  is 
shown  below. 
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vciw.py  vciin.w  vciout.w 


3.2  VCI  INPUT  FILE 

The  format  of  the  input  file  is  simply  variable/value  pairs,  one  set  per  line,  separated  by  tabs  or  spaces. 
An  example  is  shown  below. 


vciin.w 


#====== 

# 

# 

# 

# 

# 

#====== 

#  gvw 

#  m 

#  b 

#  d 

#  h 

#  he 

#  n 

#  delt 

#  CGF 

#  CEF 

#  CTF 


Description : 

Input  file  defining  variable  values  to  Calculate  One  Pass  Vehicle  Cone 
Index  (VCI)  for  Wheeled  Vehicles  using  function  calc_VCIl 


#  Gross  Vehicle  Weight  [lbs.] 

#  total  number  of  axles 

#  average  tire  section  width  (inflated,  unloaded)  [in.] 

#  average  tire  outside  diameter  (inflated,  unloaded),  [in.] 

#  average  tire  section  height  (inflated,  unloaded)  [in.] 

#  vehicle  minimum  clearance  height  [in.] 

#  average  number  of  tires  per  axle 
average  hard-surface  tire  deflection  [in.] 

#  1  if  tire  chains  are  used  or  0  if  not 

#  1  if  PWR  (power  to  weight  ratio)  <  10  [hp/ton]  or  0  if  not 

#  1  if  manual  transmission  or  0  if  automatic 


#  use  tab 
#========= 

gvw 

m 

b 

d 

h 

he 

n 

delt 

CGF 

CEF 

CTF 

serviceO 


or  space  separators 


15000 


2 

14 
40 
10 

15 
2 


278% 


3.3  VCI  OUTPUT  FILE 


The  format  of  the  output  file  is  simply  variable/value  pairs,  one  set  per  line,  separated  by  tabs.  An 
example  is  shown  below. 


vciout.w _ 

CEF  0.000000 

CGF  0.000000 

CTF  0.000000 

b  14.000000 

d  40.000000 

delt  1.000000 

gvw  15000.000000 

h  10.000000 

he  15.000000 

m  2.000000 

n  2.000000 

VCI1  28.675231 
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4.0  SUMMARY 


A  description  of  simplified  dynamic  (SSF)  and  mobility  (VCI)  metrics  relating  to  vehicle  design  parameters 
for  MDO  of  an  OCP  have  been  described.  Portable  Python  scripts,  able  to  perform  the  calculation  of  SSF 
and  VCI  on  a  wide  array  of  computing  platforms,  are  provided  in  the  Appendix.  The  use  of  ASCII  format 
(see  [13])  input/output  files  readily  accommodates  the  writing  of  parameter  data  values  and  reading  of 
output  metric  values  necessary  for  automation  of  the  simulation  process  flow  by  MDO  software. 
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APPENDIX  A.1  -  LATW.PY 


# ! /usr/bin/python 

import  sys 
import  re 
import  dynutil 


#=============================================================================== 

#  Description: 

#  latw.py  [input  file]  [output  file] 

# 

#  Reads  an  input  file  of  variable/value  pairs,  calculates  SSF,  and 

#  writes  an  output  file  of  the  variable/value  pairs,  including  SSF. 


#  Input  File  Format  Example : 

# 

# I #=========================================== 

#|#  t  #  track  width  [in.] 

#|#  h  #  CG  height  above  ground  [in.] 

# I #=========================================== 

# | #  use  tab  or  space  separators 
# I #=========================================== 

#  1 1  76 

#  I  h  36 

# 

#  Output  File  Format  Example : 

# 

# | h  36.000000 

#|t  76.000000 

# | SSF  1.055556 

# 


#  Modified: 

#  2012-04-05  Initial 


#  References : 

#  1)  see  dynutil. py 

# 


# 

# 

Variables : 

# 

fi 

input  file  name 

# 

# 

# 

fo 

output  file 

arr 

list  of  variable  and  values  per  line 

# 

vars 

dictionary  of  variable/value  pairs 

# 

vars  num 

minimum  number  of  variable/value  pairs  needed 

# 

# 

cnt 

count  of  variable/value  pairs  obtained 

if  (len (sys . argv)  <  3): 

print  "Usage:  latw.py  [input  file]  [output  file] 
sys  .exit  () 

fi  =  sys. argv [1] 
fo  =  sys. argv [2] 
arr  =  [] 
vars  =  { } 
vars_num  =  2 
cnt  =  0 


#=========================: 

#  Process  input  file 

#=========================: 


9 


fin  =  open(fi,'r') 

for  line  in  fin: 

1  =  line . strip ( M\n" ) 

arr  =  re. split (" [  \t]*",l) 

if  (arr [0] [0]  ==  "#") : 
continue 

if  (len (arr)  <  2) : 

print  "Error  in  Variable/Value  pairs . " 
sys  .exit  () 

else : 

vars[arr[0] ]  =  float (arr [1] ) 
cnt  =  cnt  +  1 

fin. close () 

if  (cnt  <  vars_num)  : 

print  "Only  %d  Variable/Value  pairs  found,  need  %d."  %  (cnt,vars_num) 
sys  .exit  () 


#==========================: 

#  Process  output  file 

#==========================: 


fon  =  open(fo,'w') 

for  key  in  sorted (vars . iterkeys ()) : 

print  »  fon,"%s\t%f"  %  (key , vars [key] ) 

v_t  =  vars [ " t" ] 
v_h  =  vars [ "h" ] 

ssf  =  dynutil . calc_SSF (v_t ,  v_h) 
print  »  fon,"%s\t%f"  %  ("SSF", ssf) 

fon. close () 

#============================================== 
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APPENDIX  A.2  -  VCIW.PY 


# ! /usr/bin/python 

import  sys 
import  re 
import  dynutil 


#=============================================================================== 

#  Description: 

#  vciw.py  [input  file]  [output  file] 

# 

#  Reads  an  input  file  of  variable/value  pairs,  calculates  VCIl,  and 

#  writes  an  output  file  of  the  variable/value  pairs,  including  VCIl. 


#  Input  File  Format  Example : 

# 


#=============================================================================== 

#  gvw  #  Gross  Vehicle  Weight  [lbs.] 

#  m  #  total  number  of  axles 

#  b  #  average  tire  section  width  (inflated,  unloaded)  [in.] 

#  d  #  average  tire  outside  diameter  (inflated,  unloaded) ,  [in. 

#  h  #  average  tire  section  height  (inflated,  unloaded)  [in.] 

#  he  #  vehicle  minimum  clearance  height  [in.] 

#  n  #  average  number  of  tires  per  axle 

#  delt  #  average  hard-surface  tire  deflection  [in.] 

#  CGF  #  1  if  tire  chains  are  used  or  0  if  not 

#  CEF  #  1  if  PWR  (power  to  weight  ratio)  <  10  [hp/ton]  or  0  if  not 

#  CTF  #  1  if  manual  transmission  or  0  if  automatic 

#=============================================================================== 

#  use  tab  or  space  separators 

#=============================================================================== 

gvw  15000 

m  2 

b  14 

d  40 

h  10 

he  15 

n  2 

delt  1 

CGF  0 

CEF  0 

CTF  0 

Output  File  Format  Example : 


#  |  CEF 

0.000000 

#  |  CGF 

0.000000 

#  |  CTF 

0.000000 

#  I  b 

14.000000 

#  |  d 

40.000000 

# 1 delt 

1.000000 

#  1  gvw 

15000.000000 

#  I  h 

10.000000 

#  |hc 

15.000000 

#  1  m 

2.000000 

#|n 

2.000000 

#|VCI1 

28.675231 

#======================================: 

# 

#  Modified: 

#  2012-04-05  Initial 


# 


# 


References : 

1) 


see  dynutil . py 


# 

#  Variables : 

#  fi  input  file  name 

#  fo  output  file 

# 

#  arr  list  of  variable  and  values  per  line 

#  vars  dictionary  of  variable/value  pairs 

#  vars_num  minimum  number  of  variable /value  pairs  needed 

#  cnt  count  of  variable/value  pairs  obtained 

# 

#===========================================================: 


if  (len (sys . argv)  <  3): 

print  "Usage:  vciw.py  [input  file]  [output  file] 
sys  .exit  () 

fi  =  sys. argv [1] 
fo  =  sys. argv [2] 
arr  =  [] 
vars  =  { } 
vars_num  =  11 
cnt  =  0 


#=========================: 

#  Process  input  file 

#=========================: 


fin  =  open(fi,'r') 

for  line  in  fin: 

1  =  line . strip ( "\n" ) 

arr  =  re. split ("[  \t]*",l) 

if  (arr [0] [0]  ==  "#") : 
continue 

if  (len (arr)  <  2) : 

print  "Error  in  Variable/Value  pairs . " 
sys  .exit  () 

else : 

vars [arr [0]]  =  float (arr [1] ) 
cnt  =  cnt  +  1 

fin. close () 

if  (cnt  <  vars_num)  : 

print  "Only  %d  Variable/Value  pairs  found,  need  %d."  %  (cnt ,  vars_num) 
sys  .exit  () 


#=============================================================================== 

#  Process  output  file 

#=============================================================================== 

fon  =  open(fo,'w') 

for  key  in  sorted (vars . iterkeys ())  : 

print  »  fon,"%s\t%f"  %  (key , vars [key] ) 

v_gvw  =  vars [ "gvw" ] 
v_m  =  vars [ "m" ] 
v_b  =  vars [ "b" ] 
v_d  =  vars [ "d" ] 
v_h  =  vars [ "h" ] 
v_hc  =  vars ["he"] 
v_n  =  vars [ "n" ] 
v_delt  =  vars["delt"] 
v_CGF  =  vars [ "CGF" ] 
v_CEF  =  vars["CEF"] 
v_CTF  =  vars["CTF"] 

veil  =  dynutil . calc_VCIl (v_gvw,  v_m,  v_b ,  v_d ,  v_h ,  v_hc,  v_n ,  v_delt,  v_CGF,  v_CEF,  v_CTF) 
print  »  fon,"%s\t%f"  %  ( "VCIl" , veil) 
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f on. close  () 


#= 
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APPENDIX  A.3  -  DYNUTIL.PY 


def  calc_VCIl (gvw,  m,  b,  d,  h,  he,  n,  delt,  CGF,  CEF,  CTF) : 

# 


#  Description: 

#  Calculate  One  Pass  Vehicle  Cone  Index  (VCI)  for  Wheeled  Vehicles 

#  (Lower  numbers  represent  better  mobility.  Less  than  20  is  good.) 

# 


#  Modified: 

#  2012-04-05  Initial 


# 

# 

# 

# 

# 

# 

# 

# 


References : 

1)  Annex  E  -  HMMWV  MECV  Performance  Specification 

2)  Clarification  of  vehicle  cone  index  w/ref  to  MMP 
Journal  of  Terramechanics  43  (2006)  85-96 

J.  D.  Priddy,  W.  E.  Willoughby 
(www. elsevier . com/locate/ j terra) 


# 

Variables : 

# 

gvw 

Gross  Vehicle  Weight  [lbs.] 

# 

m 

total  number  of  axles 

# 

b 

average  tire  section  width  (inflated,  unloaded)  [in.; 

# 

d 

average  tire  outside  diameter  (inflated,  unloaded) , 

# 

h 

average  tire  section  height  (inflated,  unloaded)  [in 

# 

he 

vehicle  minimum  clearance  height  [in.] 

# 

n 

average  number  of  tires  per  axle 

# 

delt 

average  hard-surface  tire  deflection  [in.] 

# 

CGF 

1  if  tire  chains  are  used  or  0  if  not 

# 

CEF 

1  if  PWR  (power  to  weight  ratio)  <  10  [hp/ton]  or  0  if  not 

# 

u 

CTF 

1  if  manual  transmission  or  0  if  automatic 

ff 

# 

w 

average  axle  loading  [lbs.]  =  gvw  /  m 

# 

CPF 

contact  pressure  factor 

# 

TEF 

Traction  Element  Factor 

# 

WLF 

Weight  Load  Factor 

# 

CF 

Clearance  Factor 

# 

GF 

Grouser  Factor 

# 

EF 

Engine  Factor 

# 

TF 

Transmission  Factor 

# 

WF 

Weight  Factor 

# 

DCF 

Deflection  Correction  Factor 

# 

MI 

Mobility  Index 

# 

veil 

Vehicle  Cone  Index  (One  Pass) 

# 

#= 


w  =  gvw  /  m 

#  Mobility  Index  Factors 

#  Contact  Pressure  Factor  (CPF) 
CPF  =  w  /  (0 . 5  *  n  *  d  *  b) 

#  Traction  Element  Factor  (TEF) 
TEF  =  (10.0  +  b)  /  100.0 

#  Weight  Load  Factor  (WLF) 

WLF  =  w  /  2000.0 

#  Clearance  Factor  (CF) 
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CF  =  he  /  10.0 


#  Grouser  Factor  (GF) 

GF  =  1  +  0.05  *  CGF 

#  Engine  Factor  (EF) 

EF  =  1  +  0.05  *  CEF 

#  Transmission  Factor  (TF) 

TF  =  1  +  0.05  *  CTF 

#  Weight  Factor  (WF) 
if  (w  <  2000.0) : 

CWF1  =  0.553 
CWF2  =0.0 

elif  (2000.0  <=  w  and  w  <  13500.0) : 

CWF1  =  0.033 
CWF2  =1.05 

elif  (13500.0  <=  w  and  w  <  20000.0) : 

CWF1  =  0.142 
CWF2  =  -0.42 

elif  (20000.0  <=  w  and  w  <  31500.0) : 

CWF1  =  0.278 
CWF2  =  -3.115 
elif  (31500.0  <=  w) : 

CWF1  =  0.836 
CWF2  =  -20.686 

WF  =  CWFl  *  (w  /  1000.0)  +  CWF2 

#  Deflection  Correction  Factor  (DCF) 

DCF  =  (0.15  /  (delt  /  h) )  **0.25 

#  Mobility  Index  (MI) 

MI  =  ( ( (CPF  *  WF)  /  (TEF  *  GF) )  +  WLF  -  CF)  *  EF  *  TF 

#  Vehicle  Cone  Index  (One  Pass) 
if  (MI  <=  115.0) : 

VCI1  =  (11.48  +  0.2  *  MI  -  (39.2  /  (MI  +  3.74)))  *  DCF 
elif  (MI  >  115.0) : 

VCI1  =  (4.1  *  MI  **  0.446)  *  DCF 

return  VC II 


#============== 

# 

def  calc_SSF(t, 


# 


h)  : 


#  Description: 

#  Calculate  Static  Stability  Factor  (SSF) 

#  (Higher  values  represent  better  stability.  1.0  is  good.) 


#= 

# 

# 

# 

# 

#= 


Modified: 

2012-04-05 


Initial 


#  References : 

#  1)  Special  Report  265  -  "The  National  Highway  Traffic  Safety 

#  Administration's  rating  system  for  rollover  resistance  :  an 

#  assessment",  Transportation  Research  Board  (TRB) .  ISBN  0-309-07249-2 


# 

# 

# 

# 


2)  http://www.nhtsa.gov 

3)  http: //www. safetyresearch.net/safety-issues/rollover-stability 


# 
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Variables : 

t  track  width  [in.] 

h  CG  height  above  ground  [in.] 

SSF  Static  Stability  Factor 


#  Static  Stability  Factor  (SSF) 
SSF  =  t  /  (2  *  h) 


return  SSF 


